Two - dimensional solitons in media with the stripe - shaped nonlinearity modulation 
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We introduce a model of media with the cubic attractive nonhnearity concentrated along a single 
or double stripe in the two-dimensional (2D) plane. The model can be realized in terms of nonlinear 
optics (in the spatial and temporal domains alike) and BEC. In recent works, it was concluded that 
search for stable 2D solitons in models with a spatially localized self-attractive nonlinearity is a 
challenging problem. We make use of the variational approximation (VA) and numerical methods 
to investigate conditions for the existence and stability of solitons in the present setting. The result 
crucially depends on the transverse shape of the stripe: while the rectangular profile supports stable 
2D solitons, its smooth Gaussian-shaped counterpart makes all the solitons unstable. The double 
stripe with the rectangular profile admit stable solitons of three distinct types: symmetric and 
asymmetric ones with a single peak, and double-peak symmetric solitons. The shape and stability 
of single-peak solitons of either type are accurately predicted by the VA. Collisions between stable 
solitons are briefly considered too, by means of direct simulations. Depending on the relative velocity 
we observe excitation, decay or catastrophic self focusing. 

PACS numbers: 



I. INTRODUCTION AND MODEL 

The propensity of localized patterns (solitons) to the 
catastrophic self focusing in two-dimensional (2D) media 
with the self-attractive nonlinearity of the most funda- 
mental cubic (Kerr) type is a well-known property which 
impedes the creation of 2D solitons in optics and/or Bose- 
Einstein condensates (BECs) [Ij. It was theoretically 
predicted for Kerr media [2I] and demonstrated experi- 
mentally in photorefractive crystals, whose nonlinearity 
is saturable [3] , that the stabilization of 2D solitons may 
be provided by periodic potentials, which can be created 
as optical lattices. 

A different stabilization mechanism for 2D solitons was 
theoretically elaborated in the form of nonlinear lattices 
i.e., a spatially periodic or localized modulation of 
the nonlinearity. In BEC this setting may be induced 
by dint of the Feshbach resonances controlled through 
nonuniform magnetic fields [8|. In optics, similar media 
may be structured as liquid-crystal-filled [9] or all-solid 
[10] microstructured fibers, using combinations of materi- 
als matched in their refractive index but featuring differ- 
ent Kerr coefficients. Actually, even in theoretical studies 
the stabilization of 2D solitons by the spatial modula- 
tion of the cubic nonlinearity is a difficult problem. In 
particular, it was not possible to find stable 2D solitons 
supported by spatially periodic nonlinear lattices simi- 
lar to their linear counterpart which readily stabilize the 
solitons. Quasi- ID nonlinear lattices, i.e., nonlinearity- 
modulation profiles which periodic functions of a single 
coordinate, give rise to a tiny stability area for 2D soli- 
tons, suggesting a conclusion that this mechanism is irrel- 
evant to physical applications Q (while quasi- ID linear 
lattices easily create stable solitons in a sizeable parame- 



ter region [ll|). In fact, the only nonlinear 2D structures 
that have thus far produced positive results as concerns 
the soliton stability are represented by a circle filled with 
the nonlinear material, which is embedded into a lin- 
ear or defocusing nonlinear medium and a lattice of 
such circles [6|. In particular, it was concluded that the 
nonlinearity-modulation profile with sharp edges (such 
as those of the circles) provide for an essentially better 
stabilization of 2D solitons than smooth profiles. 

The objective of the present work is to investigate pos- 
sibilities to create stable 2D solitons by the nonlinearity 
concentrated along a single or double quasi- ID stripe, 
rather than patterned as a periodic lattice. It will be 
demonstrated that these settings give rise to a finite sta- 
bility areas for several types of 2D solitons, on the con- 
trary to the negligible stabilization region supported by 
the quasi-lD periodic lattice Q. The results will be 
obtained in a semi-analytical form using an appropri- 
ate variational approximation (VA), and, in parallel, by 
means of numerical methods. 

Following Refs. [5|-[6^, the general scaled form of the 
2D equation with the stripe-shaped modulated nonlin- 
earity can be written as 

i^t = - (1/2) (V^cc. + ^yy) + g{x) IV^I' ^. (1) 

In this work, we consider three forms of the quasi- ID 
modulation of coefficient g{x) < accounting for the 
self-attractive nonlinearity: the Gaussian, 



g{x) = -exp (-x^) , 
a single rectangular box. 



-1 at \x\ < 2, 
at \x\ > 2, 



(2) 



(3) 
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and two symmetrically set boxes: 



9{x) = 







at d-1 <\x\ < d^l, 
at \x\ > 1, \x\ < d- 



1. 



(4) 



In all these cases, the coefficients may be fixed to their 
values adopted in Eqs. ([I])- (01) by means of rescalings, i.e., 
the modulation profiles (|2]) and (|3]) have no free param- 
eters, while their double-box counterpart (j4]) depends on 
a single parameter, d > 1, which measures half of the dis- 
tance between centers of the boxes, whose width is fixed 
to be 2. Notice, that if d < 1 we are dealing effectively 
with a single box. 

In the application to BEC, Eq. ^ is the scaled Gross- 
Pitaevskii equation for the mean-field wave function, 
il){x^y^t) (5(-[6]. In optics, the same model admits two 
different realizations. With t replaced by the propagation 
distance, Eq. ([1]) may be considered as the nonlinear 
Schrodinger (NLS) equation to govern the transmission 
of light beams in the bulk linear ambient medium, with 
the transverse structure induced by an embedded non- 
linear slab [Eqs. (|2]) or (|3])], or two parallel slabs, in the 
case of Eq. (j4|). On the other hand, the same equation 
([1]), with t replaced by z and y replaced by the temporal 
variable, r = t — z/c^ where c is the group velocity of the 
carrier wave, describes the transmission of spatiotempo- 
ral optical signals (alias 2D light bullets [1]) in a planar 
hnear waveguide, with an embedded nonlinear single or 
double stripe. Accordingly, stable 2D solitons reported 
in this paper are interpreted as solitary pulses of matter 
waves in the BEC, or spatial solitons in the bulk optical 
medium, or, finally, as the light bullets guided in the pla- 
nar medium. Stationary soliton solutions are sought for 
as il){x^y^t) = e~*^^0(x,y), where ja is the chemical po- 
tential (in terms of BEC), and function (l){x^y) satisfies 
equation 



/i(/)+(l/2)i 



7(x)|(/)p(/) = 0. 



(5) 



Equation ([T|) conserves the norm, N = 
/ j \il){x^y)\^dxdy^ which will play the role of an 
intrinsic parameter of soliton-solution families. In BEC, 
N is proportional to the total number of atoms, while 
in the bulk and planar optical waveguides it measures, 
respectively, the total power or energy. Equation ([1]) 
also conserves the Hamiltonian and ^/-component of the 
momentum. 



H 



P = ^jj {ra:^-^.r)dxdy. (7) 

For the application of the VA, it is relevant to mention 
that Eq. ^ may be derived from the Lagrangian, 



L = 



dxdy. 
(8) 



The rest of the paper is organized as follows. The soli- 
tons supported by the Gaussian profile (|2]) are considered 
in Section II, where it is demonstrated that they all are 
unstable. Stable solitons are found in Sections III and IV, 
in the models based on the single- and double-box modu- 
lation profiles (|3]) and dH, respectively. In these sections, 
the VA and numerical methods are used in parallel. In 
particular, three types of solitons are found in the setting 
based on the double stripe: symmetric and asymmetric 
ones with a single peak, and double-peak symmetric soli- 
tons. All the soliton species have their stability regions 
(for the double-peak soliton, it is very narrow). In Sec- 
tion V, collisions between solitons are studied by means 
of direct simulations. The paper is concluded by Section 
VI. 



II. THE STRIPE WITH THE GAUSSIAN 
PROFILE 

We start the analysis by considering 2D solitons sup- 
ported by the modulation profile (|2]). First, we apply the 
VA based on the natural ansatz of the Gaussian type too, 
with widths a and h: 

<P{x,y) = Aexp [- (1/2) {a'x' + h'^y'')] , (9) 

whose norm is N = J J (j)'^{x,y)dxdy = ttA^ / {ab) .The 
substitution of the ansatz into Hamiltonian ([6]) yields the 
following expression, where the squared amplitude, A^, 
is eliminated in favor of N: 



-4 



ttVI + 2a2 



(10) 



The variational equations for parameters a and 6, 
dH/da = dH/db = yield a system of algebraic equa- 
tions for a and b: 

7r(l + 2a2)^/2- V27V6(l + a2) = 0, (11) 
7r6 V2 (1 + 2a2) - a^TV = 0. (12) 

Equation ([12]) can be used to eliminate parameter b: 

b = a^N/ \iT^2(l + 2a?)] . (13) 
Substituting this into Eq. (fTTj) . we derive a quadratic 



dxdy, (6) equation for 



^2. 



2\2 



^(l + a^) = (^/7V)"(l + 2a^) 



A physical (positive) root of Eq. (p!4|) is 



TV 



2VA^2 - 47r2 



1 

2' 



(14) 



(15) 



which predicts that the 2D solitons exist provided that 
the norm exceeds a threshold value, A/thr = 27r. In fact, 
the same generic feature of solitons was found in other 
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models based on the spatial modulation of the nonlinear- 
ity, both two- Q-EI and one-dimensional [12]. 

Further, we insert b and a from Eqs. (fT3|) and (fT5]) into 
Eq. ([To]) to obtain the expression for the Hamiltonian 
predicted by the VA: 



equations 



H 



N 
16^2 L 



(16) 



Using the definition of the chemical potential, /i = 
dH/dN, the Vakhitov-Kolokolov (VK) criterion [Tlj, 
dji/dN < 0, which is assumed to be a necessary con- 
dition for the stability of the solitons, is then cast in the 
form of 



d^H/dN^ < 0. 



From ([T6l) we obtain 



37V 



d^H _ 1 
dN^ ~ 8^ 



IGTT^iV^ - 3^"^ - IGtt^ 

(7V2 - 4^2)3/2 



(17) 



(18) 



It is easy to check that this expression is always positive, 
hence the VA predicts that the 2D solitons supported by 
the stripe-shaped modulation of the attractive nonlinear- 
ity with the Gaussian profile cannot be stable. We also 
calculated the matrix of second order derivatives of the 
Hamiltonian with respect to a and b and check that the 
stationary point (solution of Eqs. (pTj) and (p!2]) ) is not a 
minimum. Direct simulations (not shown here) corrobo- 
rate this prediction: no stable 2D solitons can be found 
in the model based on Eqs. ([T]) and ([2]). 



III. THE MODULATION STRIPE WITH THE 
BOX-SHAPED PROFILE 

A. The variational approximation 

Proceeding to the model with the rectangular modu- 
lation profile ([3]), we again start with the VA, using two 
different ansdtze^ namely, the same Gaussian-shaped one 
as in Eq. (j9j), and also the product of hyperbolic secants: 



(x^y) = A sech(x/a)sech(?//6). 



(19) 



whose norm is N = AA^ab. An ansatz of the latter type 
was originally used for modeling multidimensional soli- 
tons in Ref. [14]. 

The substitution of ansatz (p!9|) into Lagrangian ([8j) 
yields 



L = AA^abfi 



2A2(a2 + 62) 



Sab 



-aM^sech^ i - 
9 



3sinh ( - 

a 



sinh ( - 

a 



(20) 



Deriving the Euler-Lagrange equations from this expres- 
sion, it is possible to eliminate A and 6, ending up with 



36 



N ■ 



3a2 [24 + 4a sinh(^) + a sinh(f ) 



12a cosh^(|)coth(| 



(21) 



[2 + cosh(i)] ^a {a + 24 [4sinh(i) + sinh(f )] 

(22) 

The system of Eqs. ([2T]) and ([22]) was solved numerically. 
Curves A^(/i) for the soliton families, produced by the 
ansdtze of both types, are shown in Fig. [TJ According to 
the VK criterion, portions of the soliton families depicted 
by the sohd curves may be stable. 

The Hamiltonian corresponding to ansatz (p!9|) is 



H = 



N ( 1 
6" 



1 

62 



7V2 
36a6 



2 + sech"' 



tanh 



a 
(23) 

cf . Eq. (p!6|) . Again the stability condition obtained from 
Fig. [1] using VK criterion coincide with the fact that the 
Hamiltonian (l23l) has a local minimum. 



6.50 n 



6.25- 



N 



6.00 



5.75 



-1.4 -1.2 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 



1^ 



FIG. 1: (Color online) Curves /i(A/') for the soliton family, 
as predicted by the variational approximations based on the 
sech-sech ansatz ([19]) (red) and Gaussian (|9]) (bue), in the 
model with the single-box modulation profile Q. The VK- 
stable portions of the family are shown by solid curves. Nu- 
merically found stable solitons are shown by the chain of cir- 
cles. 



B. Numerical results 

We have performed systematic numerical simulations 
of Eqs. ([1]), (j3j), with the aim to verify the predictions 
of the VA for this model. First, stationary shapes of 
2D solitons were generated by means of imaginary-time 
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simulations. Then, their stabihty was checked using the 
spht-step Fourier-transform method for simulating the 
evolution of perturbed solitons, to which a small distur- 
bance was added. The initial perturbation changed the 
total norm of the soliton by up to 0.1% (adding pertur- 
bations at this level were sufficient to clearly distinguish 
between stable and unstable solitons). 

The simulations produce stable 2D solitons in the in- 
terval of 5.645 < N < 5.85. The comparison of the nu- 
merical findings with the predictions of the VA is shown 
in the Figured) 

The fact that the replacement of the smooth Gaussian 
modulation profile by the box-shaped one agrees with 
the above-mentioned general trend discovered in other 
models dealing with the spatially modulated nonlinearity, 
viz.^ that modulation profiles with sharp edges are more 
apt to generate stable 2D solitons than their smoothly 
shaped counterparts [5|, l6|] . 

IV. THE DOUBLE-BOX MODULATION 
PROFILE 

A. The variational approximation 

To consider the model based on the dual-box modula- 
tion profile as given by Eq. (|4]), we have started with the 



VA based on the factorized-sech ansatz similar to that 
introduced above in the form of Eq. (p!9|) . but with an 
additional degree of freedom, c, which makes it possible 
to consider solitons configurations with a spontaneously 
broken symmetry: 
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FIG. 2: (Color online) An example of a stable double-peak 
symmetric soliton, with N = 6.82, /x = —0.05, and d = 1.2. 
Left: A top view shown by means of contour plots; Right: the 
section drawn through y = 0. 

The respective expression for Hamiltonian ([6j) is: 



H 



72ab 
2 + sech^ 



sech 
d^ 



d+ 1 



tanh 



1-d 



2 + sech^ 



1-d 



tanh 



1-d 



tanh 



d-1 



2 + sech^ 



c+ 1 + d 



tanh 



1 + d 



TV 



1 

(25) 



In the model with the dual-trough linear potential, for- 
mally similar to its nonlinear counterpart (j4]), the spon- 
taneous symmetry breaking of 2D solitons was stud- 
ied in detail in Ref. [15]. In a ID model with the 
nonlinearity modulation represented by two symmetric 
delta-functions, or by a superposition of two Gaussians, 
symmetry-breaking localized states were investigated in 
Ref. £]. 

The variational equations following from expression 
([25]) are cumbersome, therefore they are not explicitly 
written here. These equations were solved numerically, 
for both symmetric (c = 0) and asymmetric (c 7^ 0) con- 
figurations. The stability condition following from the 
VA was implemented as a condition that the Hamilto- 
nian must attain a local minimum at stationary points. 
The results are presented below, in comparison with re- 
spective numerical findings obtained from Eqs. ([T]) and 



B. Basic types of two-dimensional solitons. 
Numerical results. 

Numerical solutions for stationary 2D solitons were ob- 
tained by means of two numerical methods, namely, the 
imaginary-time integration (as above), and the spectral 
renormalization method [16]. In the former case, we fixed 
N and aimed to find the corresponding value of /i, while 
the latter algorithm was used to find N corresponding to 
given chemical potential /i. 

Soliton solutions of three types have been found in this 
way: (i) symmetric modes with a single peak; (ii) asym- 
metric solutions with a single peak; and (iii) symmetric 
solitons with a double peak. Solutions with two asym- 
metric peaks have not been found. Obviously, modes of 
type (iii) cannot be predicted by the variational ansatz 
([24|) (a modification of the ansatz that would incorpo- 
rate double-peak patterns leads to an extremely messy 
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algebra) . 



1. Symmetric solitons with a single peak 

In the case oi d < 1, the two boxes in modulation pro- 
file dll actually merge into one, bringing us back to profile 
(|3]) that was considered above. Accordingly, symmetric 
single-peak solitons are found in this case. At d > 1, 
there is a small region where the VA predicts stable solu- 
tions. However, the full numerical solutions do not reveal 
any stationary modes of this type. 



2. Double-peak symmetric solitons 

In the interval of 1 < d < 1.2, the spectral renor- 
malization method [l6] has produced stationary symmet- 
ric solutions with double peaks and a very shallow local 
minimum between them, see Fig. [2] (the imaginary-time 
method does not converge in this region). These soli- 
tary modes are stable only at very small values of 
when the solitons are very broad. For example, fixing 
fi = —0.05, we have found that stable double-peak sym- 
metric solutions exist in the interval of 1.1 < d < 1.2, 
where the dependence of N on {d — 1) turns out to be 
linear, see Fig. [3l It may be relevant to note that, be- 
cause the stability of the solitons was identified by means 
of direct simulations, it may happen that the stable soli- 
tons are, strictly speaking, subject to an extremely weak 
instability. Nevertheless, they are definitely stable modes 
in terms of physical applications. 
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FIG. 3: (Color online) The norm of the stable symmetric 
double-peak solitons versus the half-distance, between the 
two boxes in modulation profile (|4|, for a fixed chemical po- 
tential, /i = —0.05. 



3. Asymmetric single-peak solutions 

The system with the double-box modulation profile, 
given by Eq. (j4]) with d > 1, can support asymmetric 
solutions with a single peak. Figures |4] and [5] display the 
form of stationary wave functions in two cases, with d 
close to and far from 1, respectively. 




FIG. 4: (Color online) The same as in Fig.O but for a stable 
asymmetric soliton with N — 5.8, /x = —0.84, and d — 1.05. 
The substrate (red) line in the bottom panel represents the 
two-box profile ([4|). 
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FIG. 5: The same as in Fig. [H but for a stable asymmetric 
soliton with N = 5.8, /i = -0.962, and d = 1.8. 

The asymmetry of solutions may be naturally charac- 
terized by parameter 

J-oo I-oo dx\(l){x, y)p - dy dx\(j){x, y)\^ 

• 

(26) 

In Fig. [6]we present curves N{fi) and for the asym- 
metric solitons in the case of d = 1.2. From these plots, 
and similar ones obtained at other values of d, it is con- 
cluded that the VA based on ansatz ([24|) provides for a 
good accuracy, in the comparison with numerical results. 



C. Stability and symmetry-breaking diagrams 

The results produced by the VA and gleaned from full 
numerical solutions are collected in the form of stabil- 
ity diagrams displayed in Figs. [7] and [H It is concluded 
that the numerically found stability areas (covered by 
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FIG. 6: (Color online) Curves showing the norm and asym- 
metry parameter, N(/j) and for asymmetric single-peak 
2D solitons in the case of d = 1.2. We display the comparison 
between the results predicted by the variational approxima- 
tion (continuous/dashed curves, which depict stable/unstable 
portions of the soliton family), and their numerical counter- 
parts (chains of circles). 



FIG. 7: (Color online) The full stability diagrams for 2D soli- 
tons in the plane of {d,N). In areas a) and b), stable solu- 
tions are, respectively, symmetric and asymmetric single-peak 
solitons (black and red areas depict the stability areas as pro- 
duced, severally, by the numerical and variational methods). 
In a very thin area c) , stable symmetric double-peak solutions 
have been found, in the numerical form. 
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black oblique lines in the figures) are shifted somewhat 
downward, against the variational predictions (covered 
by red oblique lines). Around d = 1, there is a gap in 
the numerically found stability areas for the symmetric 
and asymmetric single-peak solitons. It is extremely dif- 
ficult to fill this gap using direct numerical solutions, as 
the simulations necessary to check both the existence and 
stability of the solutions turn out to be very long in this 
case. 

Further, a symmetry-breaking diagram is displayed in 
Figs. [9l and [TOl To produce this plot, we fixed the norm of 
the solitons and varied the half-distance, between the 
two wells. This implies that we keep a constant strength 
of the nonlinearity, while making the coupling between 
the two box-shaped stripes weaker. The norm of the so- 
lutions corresponding to the chain of circles in Figs.[9]and 
[TOl which were found as numerical solutions to Eqs. ([1]), 
was set as N^um = 5.8. The constant norm of the 
variational solutions presented in the same figures was 
fixed to be slightly different, NyA = 5.95, as Fig. [1] sug- 
gests a correction ratio, 6/5.85 ^ 5.95/5.8 between the 



FIG. 8: (Color online) A zoom of the stability diagram from 
Fig. close to d = 1, where the single-box modulation profile 
Q goes over into the double-box shape (|4|. 



VA-predicted solutions and their numerically found coun- 
terparts. 

The detailed picture of the symmetry-breaking tran- 
sition, displayed in Fig. [10] demonstrates that the VA 
predicts the transition of a slightly subcritical type, with 
a narrow bistability region (coexistence of stable sym- 
metric and asymmetric states) observed at 1.015 < d < 
1.034 (cf. Ref. [15], where the subcritical transition and 
bistability were demonstrated in the 2D model with the 
double- trough linear potential). However, after reaching 
the destabilization point at <i = 1.034, the line of the 
symmetric solitons does not continue to larger values of 

as one might "naively" expect. Instead, the destabi- 
lized branch of the symmetric solitons turns back, run- 
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ning parallel to the stable one. The latter feature corre- 
sponds to the existence of the pair of stable and unstable 
solution branches in the model with the single-box modu- 
lation profile, see Fig.[TJ These features of the symmetry- 
breaking diagram somewhat resembles those reported in 
Ref. ^] , which was dealing with a 2D model combining 
the double-trough linear potential and the cubic-quintic 
nonlinearity. As concerns the stable and unstable VA- 
predicted branches of the asymmetric solutions, it may 
be expected that they will meet in the limit of large 
as in that limit each of the two boxes becomes equivalent 
to the one in the model with the single-box modulation 
profile. 
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FIG. 9: (Color online) The symmetry-breaking diagram, dis- 
played as the u{d) curve for a fixed norm of symmetric and 
asymmetric single-peak solitons, see the text. The red (con- 
tinuous) and blue (dashed) curves depict, respectively, stable 
and unstable solutions predicted by the VA, while chains of 
circles represent numerically found stable solutions. 

Numerical solutions for symmetric solitons have also 
been found at some values of d > 1. However, the corre- 
sponding points are not shown in Fig. [TOl as it was too 
difficult to resolve their stability. 

Another representative stability and symmetry- 
breaking diagram can be constructed by varying the 
norm, N (i.e., the nonlinearity strength, as a matter of 
fact) at a fixed value of the half-distance d between the 
boxes. This diagram and a zoom of its most essential 
area are displayed in Figs. [TT]and[T2j The narrow inter- 
val of values of TV in which stable solitons, supported by 
the spatially modulated nonlinearity, may be found, is a 
generic feature of 2D models [4-6] . A noteworthy feature 
of this diagram is the absence of symmetric single-peak 
solutions (while their stable double-peak counterparts ex- 
ist in a very narrow interval). The stable asymmetric soli- 
tons are obviously generated by saddle-node (subcritical) 
bifurcations. Unlike the ID model with the double- stripe 
nonlinearity modulation , here the existence regions for 
the stable solitons is bounded from above, in terms of 
by the catastrophic self focusing. 



FIG. 10: (Color online) A blow-up of the symmetry-breaking 
region from Fig. [9l 
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FIG. 11: (Color online) The symmetry-breaking and stability 
diagram represented by the iy{N) curve for fixed d = 1.2. The 
red (solid) and blue (dashed) curves depict, respectively, sta- 
ble and unstable subfamilies of asymmetric single-peak soli- 
tons, as generated by the VA. The segments of black curves at 
the top and bottom represent stable solitons of the same type, 
obtained in the numerical form. The short black segment at 
the center depicts a family of stable symmetric double-peak 
solutions. 



V. COLLISIONS BETWEEN MOVING 
SOLITONS 

In the last paragraph we report our studies on the soli- 
ton collisions. In this case soliton are moving along free 
direction, transverse to the nonlinearity modulation di- 
rection. We consider two collision scenarios. One occurs 
when colliding entities are moving in single rectangular 
box, and the other when they move along two different 
boxes in a double-box modulation profile. In both case we 
observe the same diagram of the collision results. For ve- 
locities above a critical value the solitons come out from 
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FIG. 12: (Color online). A zoom of the part of Fig. II II which 
shows in detail the variational (the red curve) and numerical 
(the black curve) families of stable asymmetric single-peak 
solitons. 
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the collision being excited - the collision is inelastic. The 
process of excitation is accompanied by some radiation. 
The amplitude of the excitation grows as we approach 
critical velocity. If we go below the critical velocity the 
solitons do survive the collision, but soon after they pass 
though each other, excitation turns into instability and 
the soliton gradually dissolve. If we decrease the veloc- 
ity even further we eventually enter the region when the 
interaction time is long enough for the catastrophic self 
focusing to arise. 

First we considered the case when nonlinearity modu- 
lation function is of the form of a single rectangular box 
(see Eq. (j3j)). Stable solitons in this case exist only for 
the value of the norm between 5.645 and 5.85 as we know 
from previous considerations (see Sec. Ill, Fig.[T]). Three 
generic type of dynamics described above are shown in 
Fig. [131 The plot of critical velocity versus norm of each 
of the soliton is presented in Fig. [TU The nonmonotonic 
behavior can be understood as know that single soliton 
is least stable at the borders of its existence region. In 
the case of double-box modulation profile we in principle 
see the same kind of phenomena. We collide two asym- 
metric solitons which principal part is occupying different 
channels. The results are presented in Fig. [151 



VI. CONCLUSIONS 

In this work, we have considered 2D solitons in the 
model of the medium with the self-attractive nonlinear- 
ity whose strength is subject to the spatial modulation 
in the form of a single or double stripe. The model 
may be implemented in BEG and nonlinear optics. It 
is known that the stabilization of 2D solitons by means 
of the spatial modulation of the self-attractive nonlinear- 
ity is a challenging problem. Using the combination of 
the variational approximation and direct simulations, we 
have found that the ability of the stripe-shaped modu- 
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FIG. 13: Gollisions of solitons in a single box. Numerical data 
was obtained for N = 5.8. Sections of the wavefunction along 
y axis are shown for (a) v = 10 (excitation), (b) v = 2 (decay) 
and (c) V = 0.5 (catastrophic self focusing) respectively. 
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FIG. 14: Gritical velocity versus norm in the case of the single 
box collision. 



lation to stabilize the solitons crucially depends on its 
shape: while a smooth modulation profile in the form of 
a Gaussian does not support any stable 2D soliton, the 
rectangular (box-shaped) profile, as well as its double- 
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FIG. 15: Contour plots illustrating soliton collisions in a dou- 
ble box. Numerical results were obtained for N = 5.8 and 
d = 1.2. Three figures in the each row show contours of the 
wavefunction (a) before, (b) during and (c) after the colli- 
sion. The top row corresponds to ^; = 4. The solitons get 
excited because of the collision. The middle row corresponds 
to V = 1, where the solitons decay after the collision. Bottom 
row corresponds to v = 0.5, where we observe catastrophic 
self focusing. 



box counterpart, give rise to stable solitons. In the case 
of the double stripe, three species of stable soliton so- 
lutions were found, viz.^ the single-peak symmetric and 
asymmetric solutions, and double-peak symmetric ones. 
Single-peak solitons of both symmetric and asymmetric 
types are accurately described by the VA, and their sta- 
bility is adequately accounted for by the VK criterion. 
The symmetric double-peak solitons were found in a nu- 
merical form, turning out to be stable in a small region. 
Collisions between stable solitons were studied too, by 
means of direct simulations. We studied collisions in sin- 
gle and double rectangular boxes. Depending on the col- 
liding partners velocity we observed three types of be- 
havior: emerging solitons get some excitations for large 
velocity, for smaller velocity they decay (excitations were 
large enough to cause instability), and for even smaller 
velocity we enter the region where catastrophic self fo- 
cusing occurs. 

A challenging direction for the extension of the anal- 
ysis reported in this work is to search for stable solitons 
supported by spatial modulations of the nonlinearity in 
the 3D space. 
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